A Measurement of Time- Averaged Aerosol Optical Depth using 
Air- Showers Observed in Stereo by HiRes 
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ABSTRACT 

Air fluorescence measurements of cosmic ray energy must be corrected for attenuation of the 
atmosphere. In this paper we show that the air-showers themselves can yield a measurement of the 
aerosol attenuation in terms of optical depth, time-averaged over extended periods. Although the 
technique lacks statistical power to make the critical hourly measurements that only specialized 
active instruments can achieve, we note the technique does not depend on absolute calibration of 
the detector hardware, and requires no additional equipment beyond the fluorescence detectors 
that observe the air showers. This paper describes the technique, and presents results based on 
analysis of 1258 air-showers observed in stereo by the High Resolution Fly's Eye over a four year 
span. 

Subject headings: HiRes, extensive air-showers, atmosphere, aerosols, aerosol optical depth 
1. Introduction 
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Fluorescence detectors use the atmosphere 
calorimetrically to measure the energy deposited 
by extensive air-showers. Ultra-violet fluorescence 
emitted by particle cascades can be observed tens 
of kilometers away by a photosensitive detector 
when the primary cosmic particle is above 10 18 eV. 
The energy of the primary particle is measured in 
proportion to the total number of photons yielded 
by the shower. 

Monitoring atmospheric clarity is required to 
calibrate for atmospheric propagation losses of 
light between the shower and the detector. Ob- 
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taining this calibration requires routine measure- 
ments by specialized equipment, generally lasers, 
and LIDARS. While essential, this equipment is 
challenging to construct, maintain, and calibrate, 
especially in the remote deserts where fluorescence 
detectors are located. Active systems are lim- 
ited in their beams can not be so bright as to 
swamp the fluorescence detectors and cause satu- 
ration of the data acquisition systems. Additional 
methods to measure the aerosol optical depth and 
cross-check these conventional measurements can 
be helpful, especially when no additional equip- 
ment is needed. 

The High Resolution Fly's Eye observatory 
(HiRes), located at Dugway, Utah, USA features 
two fluorescence detector stations separated by 
12.6 km. (See Abu-Zayyad (2000) and Boyer 
(2002).) Each station views nearly the full az- 
imuth. The HiRes-1 station has one ring of tele- 
scopes that view 3.5 to 16 degrees of elevation. 
A second ring of telescopes extends the elevation 
coverage of the HiRes-2 station to 30 degrees. 
Each telescope features a 3.75 m 2 mirror that fo- 
cuses light onto a camera of 256 photomultiplicr 
tubes (PMTs). Each PMT views approximately 
(l°xl°). 

The atmosphere is modeled using molecular 
scattering and ozone absorption as a baseline. The 
remaining attenuation is attributed to aerosols. 
The HiRes experiment includes steerable lasers 
used to measure aerosol attenuation. For more 
details, see Abassi (2005). 

This paper describes an independent measure- 
ment of aerosol optical depth that uses air-showers 
viewed in stereo. It has the advantage that it is in- 
sensitive to the absolute photometric calibration of 
the detector hardware and the total fluorescence 
yield that are two of the largest uncertainties of 
the fluorescence technique. In the sense that air- 
showers are a natural part of the primary data 
sample, the technique incurs no additional cost. 
Furthermore, the measurement is made over the 
band of wavelengths that air-showers produce, and 
for the range of distances over which HiRes mea- 
sures air-showers. Nor does this analysis require 
a comprehensive reconstruction of the air-shower 
light profile, energy, or primary particle composi- 
tion; it is enough to reconstruct the shower axis, 
identify segments of the shower viewed in common 
by two detectors, and apply a set of selection cri- 



teria. 

We note that the technique has limitations. 
The relatively low flux of extensive air-showers re- 
stricts the statistical power of the technique to the 
measurement of one parameter, total aerosol op- 
tical depth, averaged over years. To reduce sensi- 
tivity of the result to details of the aerosol vertical 
distribution, the technique assumes that most of 
the aerosol is distributed below the shower seg- 
ments used in the analysis. The assumption is 
supported by an analysis of laser shots (Abassi 
2005), that found that the aerosol vertical distri- 
bution is consistent with an average scale height 
of about 1 km. For this analysis we use shower 
segments at least 1 km above the detectors. 

2. Stereo light balance method 

For stereo observations of atmospheric events to 
yield consistent results between detectors, an ac- 
curate description of the atmospheric attenuation 
is required. Conversely, a consistency constraint 
can be used to find the total optical depth. 

Here we use a molecular description of the at- 
mosphere as a baseline and apply a consistency 
constraint to find the remainder optical depth due 
to atmospheric aerosols. 

2.1. General solution 

The aerosol atmosphere is modeled with a total 
aerosol optical depth t. A ray traveling vertically 
to infinity is attenuated by one exponent of r. A 
inclined ray traveling to an altitude z is attenuated 




Detector 1 



Fig. 1. — Diagram of a cosmic ray air-shower as 
viewed in stereo. 
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(1) 

where T is transmission, a is the elevation angle of 
the ray, and SH is the Scale Height of the aerosol 
distribution. 

A useful parameter r\ can be factored out of this 
expression. 



V = 



sin a 



1 



(2) 
(3) 



Equation 3 gives us a particular definition of r\ ap- 
plicable to this model for monochromatic light in 
the aerosol atmosphere. However, the argument 
that follows requires only that rj be a known pa- 
rameter that satisfies Equation 2. 

Suppose that a segment of an extensive air- 
shower (Figure 1) produces N optical photons and 
some number, Sd, are recorded by a fluorescence 
detector during atmospheric conditions that are 
less than perfectly clear (i.e. through aerosol 
haze). Had conditions been perfectly clear (i.e. 
molecular with no aerosol) a greater fraction of the 
photons produced would have reached the detec- 
tor. Thus the same measured value of Sd would 
have corresponded to a smaller number, Nm, of 
photons produced, where Nm < N. 

Nm depends on the detected signal and, by defi- 
nition, does not depend on the aerosol property we 
wish to measure. It can be calculated using Nm = 
Sd * /, where / is a function of the measured 
shower-detector geometry, and molecular optical 
depth. The latter can be calculated from molec- 
ular scattering theory and knowledge of the at- 
mospheric density profile derived from radiosonde 
data. N and Nm are related by N — Nm/T, ig- 
noring multiple scattering effects. 

When two detectors observe the same shower 
segment, two simultaneous equations can be writ- 
ten. 



N{1) = 7V M (l)e (T " (1)) 
N(2) = N M {2)e( Tr,(2)) 



(4) 
(5) 



We constrain the two detectors to agree on the 
number of photons emitted, N(l) = N(2), and 
solve to find the light balance equation 



where 'light balance' A^v = In ( ^[2] ) and 'event 
asymmetry' A v = 77(2) — rj(l). 

It follows that a plot of A at versus A,, for a 
sample of events will have a slope of r. 

2.2. Polychromatic approximation 

Any application of the light balance equation 
(6) requires a definition of r/ which fits the atten- 
uation model involved and satisfies Equation 2. 
Equation 3 provides a definition of 77 which is suit- 
able if the light source is monochromatic. However 
for polychromatic light in the atmosphere it is not 
possible to satisfy Equation 2 with a simple defi- 
nition of r\. This difficulty is rooted in the wave- 
length dependence of scattering. Aerosol optical 
depth t is approximately inversely proportional to 
wavelength. Molecular scattering follows a much 
steeper relationship. 

For convenience, we will refer to the aerosol op- 
tical depth at 355 nm as r = t(355), with the 
understanding that depths at other wavelengths 
can be scaled from this value. 

A good approximate solution is found by mak- 
ing a guess t' close to r and redefining 



7]= - 



(7) 



tA 



(6) 



where N' is the number of photons emitted by 
the event assuming r', and Nm is the number of 
photons emitted assuming a molecular atmosphere 
(not including t'). 

This definition rigorously satisfies Equation 2 
only if t' = r, and in general it may be necessary 
to apply this solution iteratively to converge on a 
value of r, unless the approximation is particularly 
good. The quality of the approximation depends 
primarily on the spectral bandwidth. 

In the study that follows, the sensitivity of r 
with respect to r' is less than 1:100. We will set 
t' to 0.040 for the remainder of the discussion, 
since this uncertainty is much smaller than other 
errors in the analysis. 

2.3. Line sources 

Equation 6 applies to point sources of light in a 
straightforward fashion. An air-shower has a cross 
section hundreds of meters wide and is observed 
more than 10 km away, and can be considered a 
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point source traveling at the speed of light. A 
simpler approach in practice, is to treat the air- 
shower as a line source. 

A line source can be treated by integrating an 
infinite number of point sources along the line seg- 
ment. Equation 6 can be applied this way, pro- 
vided that A,, is relatively constant over the seg- 
ment. 

In this analysis, air-shower tracks are split until 
A,, varies by 0.3 or less over the track segments. 
Detector pixel size is sometimes a limitation in 
splitting the tracks. If the variation in A I( can't 
be kept below 1.0, the event is removed from the 
data. 

2.4. Data Selection 

Data from the HiRes detectors is matched by 
trigger time to produce stereo candidates. These 
candidates are passed through a Rayleigh filter 
to select track-like events while removing various 
noise triggers. Candidates may also be cut if a 
shower-detector plane can not be fit. We start 
this analysis with 5217 stereo shower candidates 
collected between December 1999 and December 
2003. 

From the stereo candidates we select 1258 
events for light balance analysis that have well 
reconstructed geometries and a common region 
observed by both detectors. Depending on the 
length of the common region, it may be divided in 
to segments. Tabulations for real and simulated 
data are provided in Table 1 and the selection cri- 
teria is described below. 

1) Sometimes the two detectors view different seg- 
ments of track. These events must be cut, since 
there is no overlapping region. 

2) Events that saturate the high gain FADC chan- 
nels at HiRes-2 are dropped. 

3) A random walk model is used to remove noise 
events. 

4) The reconstructed trajectory is required to be 
downward. 

5) If a track is very short, there is large uncer- 
tainty in the shower-detector plane and therefore 
a potentially large uncertainty in stereo geometry. 
Tracks are required to cover at least 4 degrees in 
each detector. 

6) If the opening angle is small between the two 
shower-detector planes, then there is large un- 



certainty in the intersection. Events with plane 
angles less than 8 or larger than 172 degrees are 
cut. 

7) Asymmetric Cherenkov scattering is a concern 
when the track is viewed at an oblique angle. 
Viewing angles below 30 or above 165 degrees are 
cut. 

8) To minimize the effects of any potential asym- 
metries, the maximum difference in viewing angle 
is set at 50 degrees. 

9) To reduce Cherenkov contribution and place 
observations above most aerosol, the segment alti- 
tude must be greater than 1000 m above detectors. 

10) Equation 6 can be applied to a linear track 
provided that A^ is approximately constant over 
the track. Track segments are cut if the variance 
in A,, across their length is greater than 1.0. 

3. Systematic Error 

Systematic errors in this analysis can not re- 
sult from calibration uncertainties of the detector 
hardware in the following sense. A wavelength in- 
dependent calibration scalar, k, applied to Equa- 
tion (5a). 

h i kx ^0t) = r <*-*> < s *> 

becomes an additive constant (In (k)). This error 
would shift the points in Figure 2 up or down, but 
would not alter the slope (r). k could represent 
an error in overall gain in one or both HiRes de- 
tectors, for example. A time dependent shift in 
calibration could smear the data vertically thus 
reducing the sensitivity of the slope measurement. 

Systematic errors can arise from effects that 
correlate with the aerosol optical path asymmetry 
A v . In this regard, we have examined the sensi- 
tivity of the slope measurement to a number of 
sources of uncertainty. Their sum in quadrature 
is 0.014 (see Table 2). 

The aerosol vertical distribution is modeled 
with a 1.0 km scale height, motivated by HiRes 
laser measurements Abassi (2005). A variation 
in the average scale height by ±0.3 km shifts the 
value by ±0.008. To estimate the effect of the 
Cherenkov light on the measurement, we gener- 
ated a sample of simulated showers without the 
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Table 1 

Selection of real and simulated stereo data. 





Filter 


Events 


Segments 


MC events 


MC segments 


0. 


Starting sample (sec text). 


5217 




9474 




1. 


Require commonly viewed segment (s) 


2219 


5492 


5095 


12174 


2. 


High gain channel is not saturated. 


1876 


4440 


4211 


9694 


3. 


Raylcigh filter 


1810 


4290 


4199 


9674 


4. 


Track trajectory is downward. 


1803 


4276 


4162 


9619 


5. 


Track length is 4 degrees in each detector. 


1738 


4211 


3979 


9436 


6. 


Stereo plane opening angle is 8 - 172 degrees. 


1683 


4089 


3834 


9090 


7. 


The viewing angle is between 30 and 165 degrees. 


1572 


3413 


3458 


7222 


8. 


The viewing angle asymmetry is less than 50 degrees. 


1472 


3200 


3230 


6625 


9. 


Altitude of segment is above 1000 m. 


1362 


2656 


2929 


5462 


10. Segment spans less than 1.0 in A^. 


1258 


2503 


2696 


5142 



Chercnkov component, r changed by 0.009. 

A number of other effects were also investi- 
gated. To estimate sensitivity to wavelength de- 
pendence effects in detector response, the shower 
data was reanalyzed with the calibration scaled by 
±10% per 100 nm. The difference in r was found 
to be =F0.002. The analysis and simulation used a 
fluorescence spectrum derived from the measure- 
ments of Kakimoto (1996) and the compilation of 
A. Bunncr (1967). Using the more recent spec- 
trum of Nagano (2004) shifts r by 0.002. Uncer- 
tainty in the shower axis geometry arising from 
the shower-detector plane resolution contributes 
an error of less than 0.002 to r. To estimate ef- 
fects of light transmission via atmospheric multi- 
ple scattering, the data was reanalyzed with an 
estimated contribution to each shower using the 
formalism of Roberts (2004). The shift in r was 
0.003. Finally, we include an estimated error of 
0.005 that arises from the non-linear response of 
an older model preamp used in some of the HiResl 
mirrors. 

4. Results 

Figure 2 shows plots of A^v versus A,, for real 
data and three Monte Carlo sets. Each point cor- 
responds to a segment of track. 

The data is binned in A,,. Each bin is fit to a 
Gaussian to determine a mean and a in A^r. The 
mean values each bin are weighted by the num- 
ber of entries and fit to a straight line. Statistical 
uncertainties are quoted with the slopes. 

Monte Carlo data is generated using random ge- 
ometries and primary particle energies, with lower 
energies weighted to approximate the measured 



Real Data 
Slope = 0.0424 
±0.0056 
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MCt a =0.001 <l 
Slope = 0.0079 
±0.0021 



MC t a =0.040 <3 
Slope = 0.0479 
±0.0017 



MC T A =0.083 O 
Slope = 0.0868 
±0.0016 



Fig. 2. — Plots of An (light balance) versus A,, 
(aerosol optical path asymmetry) for real data, 
and three Monte Carlo sets. According to equa- 
tion 6, the slope of each plot should be equal to r. 
The MC sets are generated using a constant r of 
0.0010, 0.0400, and 0.0833. Linear fits are made 
to a profile histogram, shown in black. 
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Table 2 
Systematic error estimates. 



Effect 


Approximate error to r 


Ccrcnkov contribution 


0.009 


Vertical aerosol distribution 


0.008 


Prcamp non-linearity 


0.005 


Multiple scattering 


0.003 


Detector wavelength dependence 


0.002 


Geometric reconstruction 


0.002 


Fluorescence spectrum 


0.002 


Quadrature sum 


0.014 



HiRcs energy spectrum. Three simulated sam- 
ples are generated using three different values of 
t. These are listed with the results in Figure 2. 
Identical reconstruction and analysis are applied 
to real and simulated data samples. 

The resulting fit for the real data yields an av- 
erage t of (0.042 ± 0.006(siai) ± 0.014(sys)) 

4.1. Cross check 

ft is simple to check the accuracy of an atmo- 
spheric model by plotting light balance as a func- 
tion of path difference (A r = r(2) — r(l)). (see 
diagram in Figure 1.) If the atmospheric model is 
accurate, the slope should be zero. Figure 3 shows 
real data reconstructed with three model atmo- 
spheres (r = 0.001, 0.040, and 0.100). A positive 
slope will indicate a deficit in the an optical depth 
of the model, while a negative slope indicates an 
excess. 

The plot with t = 0.04 has a small positive 
slope, indicating a r somewhat larger than 0.04, 
which is consistent with 0.042 shown in Figure 2. 

5. Conclusion 

Stereo cosmic ray showers are used to measure r 
in a manner that is independent of absolute detec- 
tor calibration. While it cannot replace the hourly 
and daily measurements obtained by specialized 
equipment, this method gives a cross check of the 
amount of aerosols present as averaged over ex- 
tended periods. This technique may be of use to 
other experiments that measure air-showers with 
more than one fluorescence detector station. A 
trade-off between sensitivity and statistics is ex- 
pected, depending on the distance between sta- 
tions. We note, in passing, that this work is the 



T A (model) = 0.001 
Slope = 0.0176 
±0.0027 
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Fig. 3. — Plots of A/, (light balance) versus A r 
(path difference) for real data using three different 
model atmospheres (r = 0.001, 0.040, and 0.100). 
The model with r = 0.040 results in the smallest 
slope. 
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first reported systematic use of air-showers to es- 
timate atmospheric clarity. 
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